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We report experiments on the overall phase diagram of granular flows on an incline with emphasis 
on high inclination angles where the mean layer velocity approaches the terminal velocity of a single 
particle free falling in air. The granular flow was characterized by measurements of the surface 
velocity, the average layer height, and the mean density of the layer as functions of the hopper 
opening, the plane inclination angle and the downstream distance x of the flow. At high inclination 
angles the flow does not reach an X— invariant steady state over the length of the inclined plane. For 
low volume flow rates, a transition was detected between dense and very dilute (gas) flow regimes. 
We show using a vacuum flow channel that air did not qualitatively change the phase diagram and 
did not quantitatively modify mean flow velocities of the granular layer except for small changes in 
the very dilute gas-like phase. 

PACS numbers: 45.70.Mg, 45.70.-n 



I. INTRODUCTION 



Granular flows are often separated into two categories 
dealing with (i) dense granular flows where the flow den- 
sity is not far from the density of static packing (for stud- 
ies in 2 dimensions see , H, S B or in 3 dimensions 
see d, 0, H, [IJlS El 03) and (ii) a rapid, very dilute 
(gas) regime 11311 . which is often analyzed in terms of a 
kinetic theory [141]). Although flow on a rough inclined 
plane has become a model experiment because of its sim- 
ple geometry most of the available research focuses on the 
dense regime while less data is devoted to characterizing 
the gas regime and the transition between the dense and 
very dilute phase regimes. 

A granular layer of thickness h on a rough inclined 
plane starts flowing only if the plane inclination 9 sur- 
passes a critical angle 9 C and stops flowing when 9 is 
decreased below an angle of repose 9 r (for reviews see: 
[lj|). Experiments in 3 dimensions typically fea- 
ture a storage container with an opening of height H 
that influences the volume flow rate of material down 
the plane [1, d, [HI HH EH. The properties of granu- 
lar flows in this system have been extensively studied 
when 9 is not much larger than the flow initiation con- 
ditions [ij [l(], EL ■ For 9 m 9 r and for small volume 
flow rates, intermittent flows are observed 11711 including 
avalanches [HI] and wave- like motion [Tol. I 111. [T^]. At 
somewhat larger 9, the grains flow uniformly with a sta- 
tistically steadyflow velocity that depends on the layer 
height [1, [13, [ill- For still higher 9, the flow proper- 
ties have not been well studied althou gh an interesting 
stripe state has been reported [1, [HI |20|. One might 
expect that the layer density would decrease as 9 is in- 
creased because of higher average flow velocity which cre- 



*Electronic address: btamas@szfki.hu 



ates larger shear rates and higher granular temperature. 
Because the ratio of the down-plane force to the normal 
force diverges at 9 = ir/2, one might also ask whether a 
dense phase with a well defined layer thickness continues 
to exist at large 9. 

Flowing granular layers can be described by a set of 
macroscopic variables that depend on the system con- 
trol parameters. The two control parameters for granular 
flow on a rough inclined plane are the hopper opening H 
and the inclination angle 9. The hopper opening largely 
influences the volume flow rate from the hopper whereas 
9 controls the balance of tangential and normal gravita- 
tional force on the layer. An additional parameter that 
is harder to vary systematically is the roughness of the 
inclined plane surface. For fixed surface roughness, the 
flow properties of the layer at fixed H and 9 can be char- 
acterized by the surface velocity, the height h and the av- 
erage density p. These quantities depend on the control 
parameters and on each other in a complicated manner. 
Part of our purpose in this paper is to understand these 
relationships. 

One particular issue for granular flows at high inclina- 
tion angles was suggested by idealized numerical simula- 
tions in which gravitational forcing cannot be balanced 
by energy dissipation mechanisms, and the flow is pre- 
dicted to accelerate (2l|, . Constitutive equations re- 
cently proposed for dense granular flows suggest that the 
effective friction coefficient saturates to a finite value pi 
for high shear rates 0, HH, [24| . In this case the flow is 
expected to accelerate for 9 above 9^ =atan(/i2). As we 
will show in the present work, in accordance with the 
data presented in [9[ the dense, non-accelerating regime 
can only be observed up to tan^/tan^ = 0.85. For the 
case of density we show, that p is observed to decrease 
substantially from its closed-packed value at plane in- 
clinations starting from tan#/tan# r = 1.45 (in our case 
tan(92/tan# r = 1.7). At very high plane inclinations 
where a very dilute gas phase is observed the equation 
of state is certainly much different than for the dense 
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flows and the above criterion is not relevant. For the 
case of velocity u, our characterization of a as a func- 
tion of downstream distance demonstrates that there is a 
healing length of order the size of our inclined plane for 
moderate plane inclinations that complicates a compar- 
ison with the predictions for an acceleration threshold. 
We conclude that because of the limited plane length 
our data are not complete enough to make a quantita- 
tive prediction of the angle up to which stationary flows 
exist. 

As demonstrated in other experiments, air can some- 
times have a profound effect on the observed behavior 
of the granular flow as in, for example, segre gat ion of 
vertically- vibrated granular materials [25|, [2y, (27|, dis- 
charging hourglasses [H, or impact studies [30l |. 
Thus, in order to characterize the properties of gran- 
ular flow at high inclination angles, it is important to 
understand how air interacts with the granular flow and 
to determine the level of fluidization and grain velocity 
when the role of air drag is no longer negligible. 

In this paper, we characterize the flow of relatively 
mono-dispersed sand particles on a rough inclined plane. 
The experimental apparatus, that allows for measure- 
ments of flow conditions as a function of air pressure, 
and a characterization of the average terminal velocity 
of individual grains with different mean sizes and ma- 
terial composition are presented in Sec. [TT1 In Sec. IIII( 
we describe our measurements of granular layer veloc- 
ity, height and density as a function of hopper opening 
H . We summarize our findings in a phase diagram that 
includes the transition to a gaseous phase. Finally, we 
conclude with some discussion of the implications and 
further opportunities for the system of granular flow on 
an inclined plane. 



II. EXPERIMENT 

We start by characterizing the granular material and 
describing the experimental apparatus used to make 
quantitative measurements of the flow. For determining 
the flow properties of the granular layer, we used sand 
that was sifted with 300 and 500 /im sieves to yield a 
mean diameter of d = 400 /im. We designate this distri- 
bution as having a mean of d = 400 ^m and a standard 
deviation of 50 (im. We also used finer sand, salt and 
glass beads to help calibrate the air drag effects. 

The experimental setup used for measurements of in- 
clined plane flows is shown in Fig. [TJ A glass plate with 
dimensions 230 cm x 15 cm was set inside a 274 cm long, 
20 cm diameter cast-acrylic tube. The leftmost 40 cm of 
the tube (full with sand in the image) serves as the hop- 
per. The surface of the remaining part (190 cm) of the 
glass plate was made rough by gluing one layer of grains 
onto it or by covering the plate with sandpaper that had 
a characteristic roughness of 190 /im (80 grit). The same 
flow regimes were observed for both surfaces. Because 
the surface of the sandpaper was slightly smoother, we 



observed a small (about 10%) increase of the flow velocity 
compared to the case of grains glued onto the plate. The 
tube was rotatable about the middle so that we could set 
an arbitrary inclination angle 8. The whole system could 
be pumped down to P « 0.5 mbar. 




FIG. 1: Color online. Photograph of the experimental setup. 

The flow velocity u at the surface was determined by 
analyzing high speed (8000 frames per second) record- 
ings. Space time plots were created by taking one line 
of the recordings parallel to the main flow. The Fourier 
transform of such an image (consisting of streaks as traces 
of particles in the flow) yields the average velocity of sur- 
face particles. The thickness of the flow h was monitored 
by the translation of a laser spot formed by the intersec- 
tion with the surface of a laser beam aligned at an angle 
of 20° with respect to the inclined plane. 

The flow properties were characterized by varying two 
control parameters: the plane inclination 34.1° < 6 < 
52.2° and the hopper opening 0.4 < H < 4.4 cm (i.e. 
10 < H/d < 110). Setting a constant H, we observed a 
slight increase of the hopper discharge rate when increas- 
ing #, e.g., see Fig. [2] The hopper discharge rate also de- 
pended on the presence of air in the system [28|, [29|, [3l| . 
In a typical hourglass geometry, when the hopper and 
the main chamber are separated, the flow of sand faces a 
counterflow of air as the hopper discharges and a pressure 
difference builds up. Without a connection between the 
two chambers in our experiment, the hopper discharge 
rate was significantly smaller (by about 50%) in the pres- 
ence of air when compared to the case in vacuum. By 
connecting the two chambers with a flexible tube, the 
counterflow was substantially reduced, and the discharge 
rate was about 80 — 90% of the discharge rate observed 
in vacuum (Fig. [2]). All the measurements presented here 
were done for the case of reduced counterflow. 

To set the stage for our discussions of air effects in 
inclined plane flows, we introduce some results related to 
the interaction of granular particles with air. One way 
in which air can affect a granular flow is through the 
drag force of air acting on an isolated particle. Given a 



3 



0.2 



0.15 



3 

B 

S 0.1 
2 0.05 

T3 



A vacuum 

A air 

O vacuum (with pipe) 

• air (with pipe) 



H= 2.8 cm 



A 
O 



O A 



& 



35- 



40 45 
(deg) 



50 



FIG. 2: Color online. Hopper discharge rate with H — 2.8 
cm vs. 6 at ambient pressure (filled symbols) and at P = 
0.5 mbar (open symbols). The two datasets correspond to 
the case of two separated chambers (triangles) and reduced 
counterflow by a connecting pipe (circles). 



spherical particle of radius R (diameter d) and density 
p p falling under gravity g in air with dynamic viscosity 
p and density pp, the equation of motion is 



dv 9pvc(v) 
~di + 2p p R 2 
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(1) 



where the buoyancy term is ignored since p„ >> pf 
and c(v) is a turbulent drag correction [331 c ( v ) = 
1 + 0.15(?;d/V) 2 / 3 ) with v = p/ pp. The turbulent drag 
coefficient c(v) yields a good approximation for the range 
of velocities that are of interest here. The terminal ve- 
locity vt is obtained by setting dv/dt = and solving 
numerically the equation 



0.15(— ) 2/3 v 5/3 
v 



2 P P R 2 9 
9 p 
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The air drag on the individual grains is characterized 
by the Reynolds number Re = vdjv of a particle of di- 
ameter d falling at velocity v. Taking spherical particles 
with R = 0.2 mm and with the density of sand p p = 2.5 
g/cm 3 , by solving Eq. ([2]) numerically we obtain vt = 3.3 
m/s with a corresponding Re = 90. 

The values of vt were also measured by dropping par- 
ticles from 2.7 m. All of the particles reached a terminal 
velocity where the change in their velocities was below 1% 
over a 7 cm interrogation window. The velocity distribu- 
tion reflecting the variations in particle size and shape is 
presented for 100 particles in Fig. [3] The average value 
of the terminal velocity for the 400 pm sand is vt = 2.7 
m/s. The measurements were also done for fine sand 
d = 200 ± 50 pm, salt d = 400 ± 50 pm and glass beads 
d = 510 ± 50 pm, yielding the the velocity distributions 
presented in Fig. [3] The corresponding mean values of 



vt are 1.8 m/s, 3.0 m/s and 4.0 m/s, respectively The 
calculated values of vt are indicated with a vertical line 
for each case. The calculated values of vt are close to 
the measured values with deviations between data and 
theory presumably resulting from non-spherical shapes 
and/or uncertainties in the mean particle diameter. For 
example, the deviations observed for the two sets of sand 
can be consistently explained with non-centered size dis- 
tributions. 
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FIG. 3: Color online. Distribution of experimental values 
of vt for different granular materials: fine sand (d — 200 ± 
50 pm), sand (d = 400 ±50 pm), salt (d = 400 ±50 pm), and 
glass beads (d = 510 ±50 pm). The dashed lines are Gaussian 
fits to the data, and the solid vertical lines are the calculated 
values of vt- 



Having determined that our calculation of vt agrees 
well with the measured data, we used Eq. [5] to estimate 
the air drag reduction and the resulting increase in vt 
as a function of air pressure. The dynamic viscosity p 
of air is almost constant as a function of air pressure P. 
Using an ideal gas relationship between pp and P and 
Eq. [2 one obtains vt (and Re) as a function of P. As 
seen in Fig. 01 vt is expected to increase by a factor 
of 4 when P decreases to 0.5 mbar. If air drag plays a 
significant role in determining the granular flow state, the 
substantial decrease in air drag under vacuum will reveal 
that effect. Here we ignore the expected further increase 
in the terminal velocity at very low pressures when the 
air mean free path is comparable to the particle diameter. 

Another approach to estimating the effect of air drag 
on a flowing granular layer is to assume that the gran- 



4 




0.0001 0.001 0.01 0.1 1 

P [atm] 

FIG. 4: Calculated terminal velocity of spherical particles 
with d — 400 /iia and p — 2.5 g/cm 3 vs. P. 

ular layer is similar to a fluid and that the air forms a 
boundary layer between air at large distance with zero 
velocity and the granular layer flowing with a "surface 
velocity" v. A Prandtl boundary layer description yields 
a boundary layer thickness S = 3A(vx/v) 1 ^ 2 where x is 
the downstream distance. The force per unit area Ja ex- 
erted on the flow by the entrained air is approximately 
0.66 (pf/J-v 3 /x) 1 / 2 . The relevant parameters of the air are 
IX = 1.789 ■ 1CT 5 Pa s, p F = 0.96 kg/m 3 (for local pres- 
sure P ps 0.76Patm), an d v = 0.19. Taking a distance 
x = 100 cm and v = 300 cm/sec, one obtains 6 « 0.2 
cm and Ja ~ 0.01 N/m 2 which is less than 0.1% of the 
gravitational force — a negligible effect. The surface of 
a granular flow, however, is not as sharply defined as it 
is for liquid flows. If there were substantial fluidization 
of the granular layer, the cumulative drag force experi- 
enced by the low density particles might play a role in 
the instabilities of the homogeneous flow. 

III. RESULTS AND DISCUSSION 
A. Phase diagram 

Based on qualitative observations of the flow and on 
quantitative measurements presented below, we can sep- 
arate the granular flow into regions with certain charac- 
teristic features. Such a phase diagram identifying the 
different flow regimes as a function of the two control 
parameters H and 9 is shown in Fig. [5] The ratios 
tan#/tan# r and H/d are also shown on the top and right 
axes, respectively, where 9 r = 30.5 is the independently 
measured bulk angle of repose for the granular material 
used here. The boundaries are neither precisely defined 
nor indicative of sharp transitions between different flow 
regimes. Further, this phase diagram does not capture 
the convective nature of the granular flow in that features 
or flow properties may evolve over length scales compa- 



rable to the length of the inclined plane. Nevertheless, 
identification of the general regimes are useful in setting 
the stage for more complex issues raised below. 
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FIG. 5: Color online. Phase diagram of the flow regimes in 
the phase space defined by H and 8. The dotted line divides 
the diagram into two regions: low flow rates and not too steep 
planes where a;— invariant steady flow is observed and higher 
flow rates and faster flows where the flow is still accelerating 
as measured at x = 1 m. The ratios tan 9/ tan# r and H/d 
are labeled on the top and right axes, respectively. 

At very slow flow rates and low plane inclinations the 
homogeneous flow is unstable and organizes itself into 
waves. The properties of such waves have been studied 
extensively |l0| . By increasing the flow rate and staying 
in the range of not very steep plane inclinations (up to 
about 40°), the flow becomes homogeneous. For steeper 
planes an instability occurs leading to a pattern consist- 
ing of lateral stripes [1, Qj| . A detailed characterization 
of the stripe state in our system will be presented else- 
where 120]. The dashed line divides the diagram into 
two regions: accelerating and x— invariant steady flows. 
Above the dashed line the acceleration of the flow (aver- 
aged over the range of 60 < x < 140 cm) was more than 
5% of g sm9. This boundary depends on the distance 
downstream x at which the acceleration was measured 
as discussed in more detail in Sec. IIII CI The phase dia- 
gram was found to be qualitatively the same for flow at 
ambient pressure and flow at low pressure with P = 0.5 
mbar. 



B. Flow thickness and density 

In this section, we describe our procedure for measur- 
ing layer height and mean density of the flowing layer. 
The flow thickness was measured by laser deflection. The 
plane of the images shown in Fig. [5] was parallel to the 
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FIG. 6: Image of the laser spot for (a) no flow and for flow 
at inclination angles 9 of (b) 35.0°, (c) 37.2°, (d) 40.0°, (e) 
42.6°, (f) 45.8°, (g) 48.5° and (h) 52.2°. H = 0.4 cm and 
x — 79 cm. 



inclined plane. A laser beam projected from the left hand 
side of the image at an angle of 20° with respect to the 
plane of the image produced a localized laser spot. The 
first image, Fig. [5^, shows the case of the empty chamber 
without flow, where the laser beam was reflected by the 
rough surface of the inclined plane. The other images 
(b-h) were taken in the presence of flow at a constant 
hopper opening H = 0.4 cm and downstream distance 
x = 79 cm with increasing 9 as we go from (b-h). The 
laser beam was reflected from the particles at the sur- 
face of the flowing layer. The horizontal shift of the laser 
spot with respect to image (a) measured the thickness of 
the flowing layer. This measurement is straightforward 
when the laser spot did not change its shape, i.e., the 
surface of the flow was well defined (Figs. [5J>e). The 
first sign of decreasing flow density can be seen in Fig.[6f 
where the shape of the reflected laser spot changed sig- 
nificantly. The spot spreading became more pronounced 
with increasing 9 (Figs. [Bg-h). 

To get a better measure of the flow thickness in this 
regime the light intensity / detected at a height z above 
the inclined plane was integrated over 200 equally-spaced 
images taken over a time period of 3.3 s. The I{z) curves 
are shown in Fig. [7] corresponding to the last four plane 
inclinations of Fig. [5] Because of the increasing fluctua- 
tions in the laser spot intensity for the higher 9, averaging 
over many images was necessary. 




Z (cm) 

FIG. 7: Color online. Time-averaged light intensity I(z) vs. 
z for 6 values indicated. H = 0.4 cm, x = 155 cm. Solid lines 
are curve fits of the form az b e~ cz . 

For the relatively dense regime (with a compact laser 
spot), the position of the center of mass of the laser spot 
was taken as the flow thickness h. In the very dilute 
regime we integrate the I(z) curves and define h as the 
height below which 80% of the flowing grains were de- 
tected. 




f(s) 

FIG. 8: Color online. Flow thickness h vs. t for plane in- 
clinations 6 = 42.6°, 45.8°, 48.5° and 52.2°. H = 1.6 cm, 
x — 155 cm. Flow was stopped at t m 4 s. 

The flow density was measured in the following man- 
ner. A stationary flow was established and was main- 
tained for about 4 seconds. The flow was then "frozen" 
by rapidly decreasing the plane inclination. The time 
evolution of the center of mass of the laser spot is shown 
in Fig. [5] for the same set of plane inclinations as in Fig. 
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FIG. 9: Color online. Velocity u (a-c) thickness h (d-f) and calculated normalized flow density p/ p s (g-i) as a function of x 
the distance measured from the hopper gate for three values of the hopper opening H = 0.8 cm, H = 1.6 cm and H — 2.4 cm. 
The symbols correspond to the plane inclinations of 35.0° (o), 37.2° (x), 40.0° (V), 42.6° (+), 45.8° (*), 48.5° (A) and 52.2° 



UJ The flow height fluctuated in the first 4 seconds with 
larger fluctuations in h for higher values of 0. As the flow 
stopped, the height dropped to the height correspond- 
ing to the static (nearly random closed packed) material. 
The height change between the stationary flow and the 
frozen state reflected the mean density change which be- 
came larger with increasing 0, see Fig. [51 The ratio of 
the static height and the height of the stationary flow 
gives the normalized depth-averaged density of the flow 
jo/ p s , where p s is the static density. This measurement of 
height and density averages over some distance in x be- 
cause the flow does not stop instantaneously. Assuming 
an exponential decay of the velocity and a total stopping 
time of about 0.5 sec from the data in Fig. [71 one obtains 
an averaging length of between 10 and 50 cm. Below 
we describe a more local albeit less precise measure of 
density derived from measurements of u and h. 



C. Flow characterization as a function of 

Granular material exits the hopper at low velocity with 
height H and a density p s that is close to that of a random 
close-packed state. The material accelerates, thins and 
becomes less dense as the interaction of the grains with 
the rough bottom surface partially fluidizes the granular 
state. At low inclination angles, the system reaches an 



an x— invariant steady state after some healing length £ 
where layer height h, mean velocity u, and mean density 
p do not change with downstream distance x. For 9 larger 
than about 40°, the flow is not stationary as a function 
of x and more complicated states are observed. The flow 
appears to become stationary on healing lengths of order 
the plane length L for large values of 9 as described below. 

All measurements of u and h were repeated at six lo- 
cations at a distance x measured from the hopper in the 
range of 20 cm < x < 187 cm. The flow velocity u and 
layer height h are shown as functions of x for three hop- 
per openings H = 0.8 cm, 1.6 cm and 2.4 cm in Figs. 

Ei-f. 

For higher plane inclinations the flow did not reach 
a stationary state over the plane length L. The evolu- 
tion of the granular flow can be partially understood by 
an analysis of u as a function of downstream distance x 
for different 9. The flow is driven by the gravitational 
force pgsm.9 and damped by dissipation forces, e.g., fric- 
tion and/or inelastic collisions. In the simplest approach 
we can assume that the flow would approach a termi- 
nal velocity uj like u(x) — Uf — (v,f — Uo)e~ x '^ where 
£ corresponds to the healing length of the flow. Fitting 
the data in Fig. [9ji-c to this form allows for the determi- 
nation of £ and the velocities uq and Uf for different H 
and 9. In Fig. [TUb . we show data for u and fits to the 
exponential form. The data are well fit by exponential 
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FIG. 10: Color online, (a) Fits to the u(x) data taken at H = 
2.4 cm. The symbols correspond to the plane inclinations 6 
of 35.0° (o), 37.2° (x), 40.0° (V), 42.6° (+), 45.8° (*), 48.5° 
(A) and 52.2° (•). (b) Healing length £ as a function of 6 for 
H = 1.6 cm and H = 2.4 cm. 



saturation. The resulting values of £ are shown in Fig. 
[TUb . as a function of 9 for H = 1.6 cm and 77 = 2.4 
cm. The healing lengths arc of order L/2 for 9 > 43°. 
Because the data show little variation with x for smaller 
9, our fits almost certainly overestimate £ in that region. 

From information about the surface velocity u and the 
layer thickness h, additional information can be inferred 
about the granular layer despite the inability to directly 
measure the z dependence of density and velocity. If we 
assume that the u(z) and p{z) profiles do not change their 
character along x, then by conservation of mass one has 
that the mass flux per unit channel width is F = Kpu 
where the overbar denotes a depth average. Numerical 
simulations [2l| suggest that for at least some range of 
9, p is almost constant over the depth. With that as- 
sumption, we have that F w ~phu. Thus, as a first ap- 
proximation, plotting F/uh as a function of x provides 
information about the downstream evolution of the mean 
density ~p. The degree to which this is a good approxima- 
tion depends on the details of p(z) and u(z). The slight 
difference in p between the values estimated this way and 
measured by a more precise method presented Sec. IIIIB1 
can result from a nonuniform (z dependent density) at 
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FIG. 11: Color online. Acceleration averaged over 60 < x < 
140 cm for H values of 0.4 cm (o), 0.8cm (x), 1.2 cm (A), 1.6 
cm (*), 2.0 cm (o), 2.4 cm (+) and 2.8 cm (V). 



faster flows. Because we do not know the absolute value 
for F, we normalize the resulting values of p assuming 
that the flow density for the lowest plane inclination was 
near to the density of the static packing p s near the hop- 
per. For the other values of 9 we assumed a 60% increase 
of the hopper flow rate between 9 = 35.0° — 52.2° for all 
values of H. The resulting normalized mean density p/ ' p s 
is shown in Figs. [9g-i as a function of x. One sees that 
the density drops rapidly as a result of increasing veloc- 
ity for higher plane inclinations especially for the case of 
lower incoming flow rates. A more precise characteriza- 
tion of the density change using the method described in 
Section HlIB[ made at one downstream location x — 155 
cm, will be presented below. Results of the two methods 
agree within the uncertainties inherent in each approach. 

Using the u(x) curves the acceleration of the flow was 
obtained by assuming that a = du/dt « udu/dx. The 
average value of a over the range 60 cm < x < 140 cm 
is presented as a function of 9 in Fig. [11] At higher 
plane inclinations, there is a significant acceleration for 
each value of H whereas for the lowest value of 9 the flow 
is x— invariant for all the H values we investigated. The 
value of H corresponding to a = 0.05$ sin# is indicated as 
a function of 9 in Fig. [5] with a dotted line as a boundary 
between the accelerating and x— invariant steady regimes 
measured at this distance from the hopper. 

In the following we analyze the velocity it, flow thick- 
ness h and dimensionless mean flow density ~p/p s as a 
function of 9 and H at the location x = 155 cm below 
the hopper. Parts of our data correspond to accelerat- 
ing non-stationary flow. Since we are reporting averages 
obtained near the channel center, we need to determine 
the degree to which those averages depend on the lateral 
homogeneity of the flow. Thus, we first indicate some 
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FIG. 12: Color online. Surface velocity profiles u(y) of the 
flow for plane inclinations of 6 = 34.1°, 37.2°, 40.0°, 42.6° and 
45.8°. The data were taken the location x — 155 cm below 
the hopper at H — 2.0 cm. 

Lateral velocity profiles u(y) are presented in Fig. [T2J 
In the homogeneous dense flow regime, the velocity is 
relatively constant over the center 12 cm of the channel 
corresponding to about 75% of the channel width. At the 
edges of the channel there are boundary layers that arise 
from the friction with the sidewalls. As the sidewalls of 
the channel have a smooth surface, friction is significantly 
less important there when compared to the rough plane. 
The transition to the stripe state (see curves in Fig. rvj] 
with 8 values of 42.6° and 45.8°) happens at around 8 = 
41°. 

The surface flow velocity u, the height of the flow- 
ing layer h and its standard deviation defined as a = 
VXZWSEiS dz are sri0wrL as a function of H in Fig. [13] 

as measured at x = 155 cm. 

The flow velocity is a monotonically increasing func- 
tion of H for all plane inclinations. In the fluid-like 
regime h is a monotonically increasing function of H as 
well, whereas a has a constant low value reflecting the 
constant shape of the laser spot. The sign of the transi- 
tion into the gaseous phase is the rapid increase of h and 
a with decreasing H at higher plane inclinations. 

The same set of data for u, h and a is presented in 
Fig. Q3] as a function of 8. Again u is a monotonically 
increasing function of 8 for each value of H. In the fluid- 
like regime (larger values of H) the slope of the u(8) 
curves is nearly constant. At smaller H the slope in- 
creases considerably when entering the gaseous regime. 
The flow thickness h decreases with increasing 8 in the 
fluid-like regime meaning that larger flow velocity results 
in smaller flow thickness. This effect is actually stronger 
than Fig. 114b indicates, as the incoming flow rate slightly 
increases with increasing 8 (see Fig. [5]). Thus, thinning 




FIG. 13: Color online. The flow velocity u, layer thickness h 
and its standard deviation vs. H at x = 155 cm. 



of the flow as a result of larger flow velocities is slightly 
counterbalanced by the growing value of the flow rate by 
increasing 8. 

The different flow regimes can best be characterized 
by the flow density. As we see in Fig. [15] the normalized 
mean density ~pj p s drops continuously with increasing 8. 
The density drop is more dramatic for smaller incoming 
flow rate than for thicker flows. 

The first regime corresponds to dense flows (slightly 
above the angle of repose) and is often characterized by 
the Pouliquen flow rule, where the depth averaged veloc- 
ity u is proportional to h 3 ^ 2 0, [l(| ■ This regime is re- 
ported to exist for the plane inclinations tan#/tan# r < 
1.45, and the dense flow can be unstable with respect 
to the formation of waves [lfj. According to our mea- 
surements the mean density in this regime is slightly de- 
creasing with 8 but always stays larger than ~pj p s = 0.8. 
According to a recent theory by Jenkins [331 ] the den- 
sity should decrease as p/p s = 1 — B-tan 6 0. A reason- 
able fit is obtained by this formula for a range of H 
up to the plane inclination tan0/tan# r < 1.5 yielding 
B = 0.52. In the second regime falling in the range of 
1.45 < tan^/tan^r < 2.4, a stripe pattern can be ob- 
served with an average density of 0.2 < ~p/p s < 0.8. 
The detailed characterization of the stripe structure is 
beyond the scope of this paper. In the third regime, the 
stripe structure disappears as the flow gets very fluidized 
with an average density of p/ p s < 0.2. A dramatic den- 
sity decrease was observed for lower incoming flow rates, 
typically at tan#/ tan# r > 1.6, yielding a gas-like phase 
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FIG. 14: Color online. The flow velocity u, flow thickness h 
and its standard deviation oh as a function of 9 at x — 155 
cm. 
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FIG. 15: Color online. Normalized mean density p/ ' p s as a 
function of for H = 1.2 cm (x), 1.6 cm (•), 2.0 cm (A) 
and 2.4 cm (o) at x = 155 cm. The dashed lines are guides 
to the eye, while the continous line corresponds to p/p s = 
1 - O.52-tan 6 0. 



where the density was less than 5% of p s . We character- 
ize this transition in more detail below. 

The surface flow velocity u as a function of height h is 
presented in Fig. [T6a for the fluid-like regime. Note that 
only the first curve (at 9 = 35.0°) corresponds to sta- 
tionary flows, the other curves (partly) fall already in the 
range of accelerating flows. Two sets of data were taken: 



(i) in the presence of air (filled symbols) and (ii) at low 
pressure P = 0.5 mbar (open symbols). In this regime 
the two sets of data match implying that air drag does 
not become important even at the fastest flows where 
the grain velocity is close to the value of the terminal 
velocity measured in free fall. In the gaseous regime, 
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FIG. 16: Color online. The flow velocity u vs. h for a set 
of plane inclinations in the presence of air (filled symbols) 
and in vacuum (open symbols) in the (a) fluid-like phase, (b) 
gaseous phase. The transition between the above regimes is 
indicated in (c) by horizontal arrows. Data taken at x = 155 
cm. The horizontal dashed lines correspond to vt- 

however, (see Fig. [TBb) particle velocities measured at 
ambient pressure are slightly smaller than at low pres- 
sure. Thus, in the gas phase the average interparticle 
distance is considerably increased and the contribution 
of air drag to the dissipation becomes important. 

The changing nature of the flow is visualized in Fig. 
116b where all the data points are presented for two plane 
inclinations for the two regimes and the transition be- 
tween them. The sharply different tendency of the u(h) 
curves in Fig. HSb and Fig. [TBb reflects the nature of 
the transition between the fluid- like and gaseous phases. 
For 8 = 48.5° the larger flow rates give rise to dense 
flows (Fig. [TBb ) where increasing H leads to increasing 
h and faster flow. At lower values of the hopper flow 
rate decreasing H leads to increasing levels of fluidiza- 
tion and to a transition to the gaseous phase. In this 
regime, see Fig. [16b . the density of the flow rapidly de- 
creases with decreasing H and the measured thickness 
rapidly increases. Contrary to trends for the fluid-like 
phase, the leftmost data points of the curves in Fig. [TBb 
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correspond to the highest hopper flow rates and the right- 
most ones to the lowest hopper flow rate. Thus, the flow 
is not simply determined by (0, h, x) as is the case for 
dense flows. For a given h, 9 and x, two solutions exists, 
one corresponding to a dense phase and the other to a 
very dilute one, depending on the flow rate. 

A recent linear stability analysis indicated, that the 
air drag could play a role in the development of stripe 
patterns in the fluid- like regime [HI, According to 
our findings the density of the flow for the stripe state 
was in the range of 0.2 <~p/p s < 0.8, corresponding to 
the regime where the role of the air drag is minor. The 
flow properties are not affected even at flow velocities 
near the terminal velocity in free fall. At the transition 
to the gaseous state the stripes disappear and the rapid 
decrease of the flow density leads to visible effects of air 
drag. 

In summary, we have presented a detailed description 
of granular flow on a rough inclined plane - one of the 
most commonly used model systems for granular dynam- 
ics - concentrating on the fast flow regime. We developed 



a method to measure the depth averaged normalized flow 
density ~pj p a and characterized the flow regimes as a func- 
tion of ~p/p s - We have characterized the transition to a 
very dilute gaseous phase that takes place by increasing 
plane inclination and decreasing incoming flow rate. By 
measuring the flow properties at ambient pressure and 
at low pressure (P « 0.5 mbar) in a vacuum flow chan- 
nel we have shown that the dissipation by the air drag 
(with respect to the other dissipational processes) is non- 
negligible only in the dilute gas-like phase. Other impli- 
cations of this work involves the possibility of linking 
changes in the flow structure (such as the stripe state) as 
a function of plane inclination to the change in the depth 
averaged normalized flow density p/p s . 
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of Sciences and the Hungarian Scientific Research Fund 
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